Theory of quantum fluctuations of optical dissipative structures and its application to 
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We present a method for the study of quantum fluctuations of dissipative structures forming in 
nonhnear optical cavities, which we illustrate in the case of a degenerate, type I optical parametric 
oscillator. The method consists in (i) taking into account explicitly, through a collective variable 
description, the drift of the dissipative structure caused by the quantum noise, and (ii) expanding 
the remaining -internal- fluctuations in the biorthonormal basis associated to the linear operator 
governing the evolution of fluctuations in the linearized Langevin equations. We obtain general 
expressions for the squeezing and intensity fluctuations spectra. Then we theoretically study the 
squeezing properties of a special dissipative structure, namely, the bright cavity soliton. After 
reviewing our previous result that in the linear approximation there is a perfectly squeezed mode 
irrespectively of the values of the system parameters, we consider squeezing at the bifurcation points, 
and the squeezing detection with a plane-wave local oscillator field, taking also into account the 
effect of the detector size on the level of detectable squeezing. 

PACS numbers; 42.50.L; 42.65.Sf 
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I. INTRODUCTION 

Quantum fluctuations in the light field observables 
have been a subject of intensive and maintained research 
since the appearance of quantum optics in the early six- 
ties, especially since the discovery of the nonclassical 
states of light and the generation of squeezed light in the 
mid eighties [T]. Quantum fluctuations are unavoidable 
as their origin lies in the impossibility of determining 
two canonically conjugate observables with a precision 
better than that allowed by the Heisenberg uncertainty 
relations. In the case of the radiation field, quantum 
fluctuations manifest in quantities such as the photon 
number, the field phase, the Stokes parameters, or the 
field quadratures 0, H, |j, H, @| . 

Although its origin could be traced back till the early 
days of quantum mechanics Q, the modern squeezing 
research goes back to the mid seventies, its foundational 
period being closed with the experimental generation of 
squeezed light in the mid eighties In a single-mode 
state the two field quadratures (equivalent to the position 
and momentum operators of a harmonic oscillator) con- 
stitute a canonically conjugate pair, the product of their 
uncertainties being consequently limited by the Heisen- 
berg inequality. When the field mode is in a coherent 
state, the uncertainty is the same for any field quadrature 
and equals that of the vacuum state. Then, a squeezed 
state is that for which the uncertainty in a particular field 
quadrature is smaller than that of the quantum vacuum, 
a reduction achieved at the expense of an increase in 
the uncertainty of the complementary field quadrature. 
These states can be generated in nonlinear processes 
(four-wave mixing, parametric down-conversion...), and 
the field quadratures can be detected in a homodyne de- 
tection experiment in which the quantum state is mixed 
with a classical (intense, coherent) local oscillator field 



i- 

Single-mode squeezing has been intensively and exten- 
sively investigated for almost three decades ^Ij]. A well 
established result is that a high degree of squeezing is ob- 
tained in nonlinear cavities near the bifurcation points, 
and that the squeezinglevel degrades as the system is 
brought far from them Q . Multimode squeezing has also 
been considered in the past. One of the most immediate 
cases is that of nonlinear cavities working in several longi- 
tudinal modes, a type of system already exhibiting noise 
reduction features specific to multimode systems. For ex- 
ample, the quantum noise suppression on the difference 
of intensities (amplitude-squeezing) of a two-mode opti- 
cal parametric oscillator above threshold [2, S} , which is 
associated to the existence of a continuous diffusion of 
the phase difference between the two modes. It is to be 
remarked that squeezing appears in this case associated 
to a collective variable, namely the intensity difference, 
and that, remarkably, the noise reduction level is inde- 
pendent of the system proximity to a bifurcation point. 

A particularly interesting case of multimode squeez- 
ing is that of solitons. Squeezing in optical fibre solitons 
{temporal solitons) has attracted much interest since the 
late eighties [sj, [lOj and is by now quite well understood, 
see e.g. [ll| for a recent review. In this problem the rel- 
evance of collective variables, such as the position or the 
momentum of the soliton, is very clear in the sense that 
these are the observables in which the behavior of quan- 
tum fluctuations is easily detectable. As we comment be- 
low, this problem has some similarities, but also strong 
differences, with the problem of cavity solitons that we 
treat here. A related problem, that of the squeezing of 
spatial solitons in Kerr media, has also been considered 
recently [S [H, H, [H . 

For our purposes, a most exciting connection is that 
existing between pattern formation in nonlinear optical 
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cavities and squeezing [l^ [T3| ■ In these systems the con- 
cepts used in the analysis of quantum fluctuations, which 
has started relatively recently, must be generalized to 
cover correlations at different spatial points p^ . 

Extended nonlinear optical systems, specially nonlin- 
ear optical cavities with large Fresnel numbers, are sys- 
tems that spontaneously display dissipative structures, 
which are extended patterns that form in the plane or- 
thogonal to the light propagation direction, through an 
spontaneous symmetry breaking. One of the new con- 
cepts that appeared when the analysis of quantum fluc- 
tuations in these systems was addressed is the quantum 
image [l^, which can be described as a precursor ap- 
pearing below threshold of the pattern that the system 
would display above threshold. The quantum image is 
not detectable at low observation frequencies as in this 
case the fast dynamics of quantum fluctuations is washed 
out. We refer the reader to existing reviews for a re- 
sume of the main researches in the field [l^ • For our 
purposes, which concern the analysis of the squeezing in 
dissipative structures, the work carried out by Gatti and 
Lugiato [20, 21] on the squeezing of extended degenerate 
optical parametric oscillators (DOPO) below threshold, 
is particularly close (see also [23 ) , as well as the analysis 
of the role played by the Goldstone mode in [2^ [IJ] . 

In this article we investigate quantum fluctuations of 
dissipative structures in the DOPO. In a first part we 
develop the theory for the calculation of linearized quan- 
tum fluctuations, in particular for the calculation of the 
squeezing and intensity fluctuations spectra. The the- 
ory is general in the sense that no particular dissipative 
structure is assumed. Moreover, although the model for 
a DOPO in the large pump detuning is used, the derived 
expressions are easily generalizable to cover other non- 
linear cavity models. Then, in a second part we apply 
this theory to the study of the squeezing properties of a 
special dissipative structure appearing above threshold, 
the so-called bright cavity soliton. 

In [2^ we already advanced a particular property of 
quantum fluctuations which is specific to pattern forma- 
tion (i.e., that is absent when the emission is homoge- 
neous in space), which can be put in short as follows: 
Due to the translational invariance of the problem (the 
position of the pattern in the transverse plane is not 
fixed when the pumping is spatially homogeneous), there 
is a particular transverse mode that is free from quan- 
tum fluctuations (in the linear approximation), a per- 
fectly squeezed mode. And this occurs irrespectively of 
the value of the parameters of the system. The only con- 
ditions are that (i) the system be translational invariant, 
and (ii) that the output field displays a pattern. This 
particular mode corresponds to the transverse linear mo- 
mentum of the pattern. 

Here we shall go beyond this result by studying the 
squeezing properties of a special pattern, the bright cav- 
ity soliton. For that we shall consider the DOPO model 
in the large pump detuning limit, as the problem is much 
simpler in these conditions because there exists an ex- 



plicit analytical solution for the localized structure. As 
stated, the article consists of two parts. In the first part 
(covering Sections II to V) we derive general expressions 
for the study of quantum fluctuations (e.g., the linear 
squeezing spectrum) without specializing to any partic- 
ular pattern. Then, the second part (Section VI) is de- 
voted to the bright cavity soliton. Finally, in Section VII 
we give our main conclusions. 



II. LINEAR THEORY OF QUANTUM 
FLUCTUATIONS OF OPTICAL DISSIPATIVE 
STRUCTURES 

In this section we present the general method that will 
allow us determining quantities related with the quan- 
tum fluctuations of optical dissipative structures (e.g. 
the squeezing spectrum) in the linear approximation. An 
outstanding feature of the method is that it circumvents 
the numerical integration of the dynamical, stochastic 
equations, which is always a problematic (and terribly 
time consuming) task. Instead, the method exploits the 
diagonalization of the linear problem, which allows sim- 
plest formal solutions. 

We shall use the DOPO with plane cavity mirrors as 
a model for presenting the method. Any other nonlin- 
ear optical cavity with plane mirrors sustaining dissipa- 
tive structures can be studied along similar ways after 
straightforward particularization of the expressions given 
below. 

We assume that the dynamical equations of the stud- 
ied system have been cast in the form of classical-looking 
Langevin equations corresponding to some coherent state 
representation. In particular we assume that a general- 
ized P representation is being used (Sec. II. A) as 
its normal ordering equivalence allows a direct computa- 
tion of measurable quantities corresponding to the fields 
leaving the cavity. Furthermore we assume that those 
Langevin equations have been linearized around a classi- 
cal dissipative structure (Sec. II. B). The method consists 
in separating the fluctuations in two classes: (i) Those 
coming from the drift of the dissipative structure, as it 
can move freely across the transverse section owed to the 
spatial translation invariance of the model (Sec. II. C); 
and (ii) Formally solving the equations corresponding to 
the remaining (" internal" ) fluctuations making use of a 
special basis (Sec. II.D), what allows deriving a general 
expression for the linearized squeezing spectrum (Sec. 
Ill), as well as for the spectrum of intensity fluctuations 
(Sec. IV). This section ends with a general result con- 
cerning the squeezing of dissipative structures (Sec. V). 



A. Langevin equations for a planar DOPO in the 
generalized P representation 

We consider the model for a type I DOPO with plane 
cavity mirrors of [26] pumped by a plane wave coherent 
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field of frequency 2u!s and amplitude fin- An intracavity 
X^^^ nonlinear crystal converts pump photons into sub- 
harmonic photons, at frequency Wg, and vice versa. Only 
two longitudinal cavity modes, of frequencies mo (pump 
mode) and tJi (signal mode), which are the closest to 2ujs 
and LOs, respectively, are assumed to be relevant. The 
cavity is assumed to be single-ended, i.e., losses occur at 
a single cavity mirror, where the intracavity modes are 
damped at rates 7„, n = 0, 1. The above frequencies de- 
fine dimensionless pump and signal detuning parameters 
through Aq = (wo - 2ws) /70 and Ai = (wi - Wg) /71, 
respectively. 

We denote the intracavity field envelope operators for 
pump and signal modes, which propagate along the z 
direction, by (r, i) and yii(r, i), respectively, where 
r — (x, y) is the transverse position vector, obeying stan- 
dard equal-time commutation relations 



(1) 



As commented we shall use a coherent state represen- 
tation in order to handle the problem. These represen- 
tations set a correspondence between the quantum op- 
erators Am{r,t) and Aj„(r, i) and the c- number fields 
Am(Y,t) and y^^j(r, i), respectively, which are indepen- 
dent but in their (stochastic) averages, which verify 
(yl+(r,<)) = (^m(r, i))*. The physical meaning of the 
stochastic average of any function of the c-number fields 
depends on the representation. In particular we shall use 
the generalized P representation [27], generalized to in- 
clude the spatial nature of the multimode problem here 
considered [1^, in which the stochastic averages corre- 
spond to quantum expectation values computed in nor- 
mal and time ordering. 

As shown in Appendix A, in the large pump detun- 
ing limit (|Ao| > 1, |Ai| ,7o/7i) the DOPO dynamical 
(Langevin) equations in the generalized P representation 
can be written as 



dt 



Ai{Y,t) = 7i [LiAi + iiAl + i^AlAl 

\ hi 

^tT^M+^^^^Ki-, t) , (2a) 
^+(r, t) = 7i (^LIA^ + ^lAi - i^yl+Mi) + 

(2b) 



d_ 
di 



where Li = -(1 + iAi) + ill^^, h = c/VSoTrh is the 
diffraction length for the signal field, — d^/dx^ + 
d^/dy^ is the transverse Laplacian operator, and we have 
introduced the real and dimensionless parametric pump 
parameter /i 



where g is the (real) nonlinear coupling coefficient, Eq. 
()79p in Appendix A, the normalized nonlinear coupling 
coefficient 
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27o7i |Ao| 



(4) 



and cr — signAo = ±1. Finally, ?7(r,t) and 77+(r,t) are 
independent, real white Gaussian noises of zero average 
and correlations given by Eqs. (j85p in Appendix A. 

Equations ([2]) are the model we shall consider along 
this paper. 



B. Dynamics of quantum fluctuations: Linearized 
Langevin equations around the classical dissipative 
structures of the DOPO 



In the classical limit (A^ 
noises being ignored), Eqs 



7o7i |Ao| 



(3) 



being interpreted as A*, 
((57)1 have the form of 
a parametrically driven nonlinear Schrodin ger equation 
(PDNLSE), first derived for the DOPO in [13, [Hi, see 
Eq. ([88]) in Appendix B. In our case a accounts for the 
cases of self-focusing {a = +1) or defocusing (cr = —1) 
of the PDNLSE, which determine the kind of dissipa- 
tive structures (DS) supported by the DOPO. These DS 
are patterns that appear in the transverse plane with re- 
spect to the direction of light propagation. We denote 
these (steady) structures by Ai (r) = ^1 (r — ri), where 
Ti = {xi , j/i) is arbitrary due to the translation invariance 
of the problem. They can be, e. g., periodic patterns or 
localized structures 30,,, 133, ISJ, [Sg] . Although in the 
second part of this paper we shall concentrate on a par- 
ticular type of DS, namely the bright cavity soliton, we 
stress here that the treatment we present below is com- 
pletely general and covers the description of quantum 
fluctuations of any stationary DS. 

The dynamics of the quantum fluctuations around any 
DS is studied by setting 

Ai (r, t) = ^1 (r - n (t)) + ai (r - n (t) , t) , (5a) 
At (r, t) ^AUr- ri (t)) + a+ (r - n (t) , t) , (5b) 

where Ai and AI are the classical stationary mean val- 
ues of the field corresponding to a particular DS (i.e., the 
stationary solutions of Eqs. ([2]) when the noise terms are 
neglected), and ai and af are the c-number fields ac- 
counting for the quantum fluctuations. Notice that the 
position of the classical solution, ri {t) ~ {xi,yi), is let 
to vary in time in order to properly describe the diffu- 
sive movement of the DS, which is excited by (quantum) 
noise. 

Linearizing the Langevin equations around the classi- 
cal solution we get, to first order in the fluctuations, the 
linearized equation of motion for the quantum fiuctua- 
tions, that read 

- « (g,^ + Gy^) + = 7i/:ai + VTih, (6) 
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where 



G 



dxiy) 



SLi is the quantum fluctuations vector 

ai(r,t) 
h is the noise vector 

h(r,<) = 

where 



ai(r, t), 
a+(r, t) 



(a^ 77(r,t) 
a* r]+{r,t) 



(7) 



(8) 



(9) 



ao = fj. + i(TK~'^Al. (10) 
Finally, the linear operator £ and its adjoint read 



C = 



al CI 



al Ci 
2icrK^^ \Ai{r)\ 



(11a) 
(lib) 



where we note that a typo has been corrected in Eq. 
(IHbl) with respect to the corresponding expression in 

In Eq. ([6|) two terms {dxaidxi/dt and dyaidyi/dt) 
have been neglected as they are of second order in the 
fluctuations. Notice flnally that Eq. ^ has the standard 
form of a set of linearized Langevin equations but for the 
flrst term appearing on the left hand side, which describes 
possible displacements of the DS on the transverse plane. 

All the information about quantum fluctuations in the 
linear approximation is contained in Eq. ([6]) . Our goal is 
then finding an efficient method for solving it at the time 
that relevant information is extracted in a transparent 
way. This is accomplished by using the eigensystems of 
the linear operators 



C. Diagonalization of the linear problem. The role 
of the Goldstone modes: Drift of the dissipative 
structure 

Our main purpose in this work is to study the prop- 
erties of quantum fluctuations around the semiclassical 
mean value corresponding to a DS. In this section we 
solve Eqs. which will allow to characterize the quan- 
tum fluctuations, in particular through the squeezing 
spectrum. With this aim, we develop a general method 
suitable for obtaining the formal solution of Eqs. ([6]), 
suitable for any system and any classical stationary DS. 

Let us assume without proof that the set of eigenvec- 
tors of the linear operators C and C\ Eq. (fTT|) . form 
a biorthonormal basis. (In the second part of this arti- 
cle we show that this is indeed the case for the bright 
cavity soliton solution, unlike the problem of conserva- 
tive temporal solitons where the set of eigenvectors must 



be completed in order to form a proper basis [HI). The 
method used to solve Eq. © consists in expanding the 
quantum fluctuations in this biorthonormal basis. 
We denote the eigensystems of C and by 

£v,(r) = A,v,(r), v,(r) = (^^^^^^^^ , (12a) 

^^w.(r) = A*w.(r), w.(r) = (^^^^^^^^^ . (12b) 

In the above and in the following expressions, the index 
i represents both the discrete and the continuous spectra 
as we do not want to overburden the notation. Note also 
that, in the following, Kronecker deltas should be un- 
derstood as suitable Dirac deltas as well as sums should 
be understood as suitable integrals when referring to the 
continuous part of the spectra. 

We define scalar product as usual 



(u|s) = J '^^^ "^(i") ■ ^(^ 



so that relation 



(wi|/:s) = A,; (Wi|s) 



(13) 



(14) 



holds for any s. Finally, all eigenvectors are assumed to 
be normalized as 



(15) 



The spectra must be computed numerically in general. 
However, it will be convenient for our purposes to state 
two general properties of the discrete spectra. These are: 

• Property 1: For any parameter set there exist Gold- 
stone modes, Vi^(iy) = Gx{y), with G^(y) given by 
Eq. (O, satisfying 



(16) 



and the associated adjoint eigenvectors, denoted by 
'Wix(iy), which verify 



^^^lx{ly) = 0. 



(17) 



This property is a consequence of the translational 
invariance of the problem, as any DS can be located 
at any position on the transverse plane. 

• Property 2: For any parameter set there exist 
eigenvectors of the adjoint problem, which we de- 
note as Vf2xi2y), verifying 



^^'^2x(2y) — -2w2a;(2jy). 

These eigenvectors are 



W2k(2i/)(i") = dx(y) 



iAi{r) 
-lAUr) 



(18) 



(19) 



as is straightforward to be checked. This prop erty 
will be associated to a perfectly squeezed mode [23 , 
as we show below. 
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Now, the linearized Langevin equation ([6]) can be 
solved by expanding the quantum fluctuations on the ba- 
sis {v^}, 



ai(r,t) ^c,(t)v,(r), 



(20) 



where the Goldstone modes have been excluded from this 
expansion as any contribution of them to ai(r, i) would 
entail a shift of the solution, which is already accounted 
for by the (still undetermined) location of the DS, by 
definition. 

First we project Eq. ^ onto wix and "Wiy, obtaining 



Xi 



yi = 



K 

7i 



;wi,ih), 

;wiy|h) , 



(21a) 
(21b) 



We see that the DS diffuses driven by quantum fluc- 
tuations as anticipated. The diffusive drift of the DS 
through a time id can be evaluated by considering the 
mean squared deviation of the position of the DS along 
this time. 



p(t) =rl(^)-rl(^-^d)• 



(22) 



The variance <^p^(i)) can be calculated from Eas. ([2T|) . 
which give the time evolution for ri = [xi, j/i), obtaining 



Xl{t) 



/7i 



dt'(wi,(r)|h(r,0>, (23) 



and analogous expression holds for yi{t), when wi^; is 
replaced by wiy in (j23p . Substituting this solution into 
Eq. (p2| , we reach the expression of the variance of p{t) , 
which is linear in time *24!|, and reads 



where the diffusion coefficient is given by 

D = 2^^Rej dV [wUv) + wlyiv)] a*{r). 



(24) 



(25) 



The knowledge of the variance allows us to evalu- 
ate the possible effects of the DS movement on the noise 
detected, e.g., in a homodyning experiment. This could 
be quantitatively important as a possible noise reduction 
could be blurred by the diffusion of the structure [25j . 



D. Formal solution to the linearized Langevin 
equations 



Now by substituting ((20)) into ([6]) and projecting onto 
{wi}, we obtain the evolution equation for the expansion 
coefficients in Eq. (|20p : 



Ci = 7i AiCi -t- {wi\h) 



(26) 



where the index i does not include the Goldstone modes, 
as commented. We write down Eq. (pS)) in the temporal- 
frequencies space. By defining the Fourier transforms 



Ztt J 



(27a) 
(27b) 

(28a) 
(28b) 



we find a simple expression in the temporal spectral do- 
main for Eg. ipS)) . that reads 



h(r,t) = i- J diu e-*h(r,c.), 
and the corresponding inverse transforms 

h{r, uj) ^ I dt e-"^*h(r,t). 



iujCi{uj) = 7iAjCi(a;) -|- VTi (r) |h(r,w)) , (29) 
which gives 

/77(w, (r) |h(r,u;)' 



Ci(w) 



(30) 



By using solution ([30| in Eqs. ([27)) . we retrieve the 
expansion coefficients Ci{t) 



C^{t) 



27r 



^1 ' duj e''^* 



w,j (r) |h(r,u;) 



(31) 



which is the formal solution of the linearized Langevin 
equations, Eqs. and . 



E. Modal correlation spectrum 

Once the time-dependent expansion coefficients Ci{t) 
are known, Eq. (PTjl . we can calculate the two-time cor- 
relations between the different modes. The knowledge 
of these correlations, and specifically, of their spectra, 
is necessary in order to characterize some properties of 
quantum fluctuations, such as the spectrum of squeezing 
or the spectrum of intensity fluctuations of the quantum 
field exiting the nonlinear cavity, as will be shown below. 

We define the correlation spectrum between two modes 
labeled by indices i and j as usual 



S,,{lu) = J dT e-'"^ {c,{t + r), c,(r)) (32) 
where the correlation 

{c4t + T),Cj{T)) = 

(c,(i + r)c,(T))-(c,(t))(c,(t)), (33) 
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and (C) is the (stochastic) average value of O. By using 
Eqs. and Sij{uj) can be written as 

w,|h(r,^)\,(w,|h(r,^)\\. (34) 



is then measured, and it turns out to be given by [2l| 



Taking into account the properties of noise, Eqs. ([SS]) . 
one obtains straightforwardly 



Wi|h(r,L^)),(wj|h(r',c^') 



27rS {lo + uj') / d^r d,j (r, r) , (35) 



where 



dtj (r, r) 

< (r) u;* (r) ao (r) + [w+ (r)] * [w+ (r)] * (r) , (36) 

so that the modal correlation spectrum Sij{uj) can be 
written as 



5,y(w) 



(A, -iw/7i) (Aj +iw/7i)' 



(37) 



where we have introduced the matrix D^j, which we call 
modal diffusion matrix (because of the similarity of Eq. 
(157)1 with an spectral matrix [1^), defined by 



D,j = J (fr d,j (r, r) , 



(38) 



Note that all modal correlations can be obtained just 
by diagonalizing the linear problem. 



III. SQUEEZING SPECTRUM IN THE LINEAR 
APPROXIMATION 



We consider now the squeezing properties of the clas- 
sical DS Ai as measured via a balanced homodyne de- 
tection experiment [21], which allows the direct measure- 
ment of quadrature squeezing. The quantum field outgo- 
ing the nonlinear cavity, to be denoted by Ai^out (r, t), is 
combined at a beam splitter with a local oscillator field 
(LOF). This LOF lies in a classical multimode coherent 
state aL (r — tl (t)) of intensity much larger than that 
of ^i,out (r, t). (The shift tl {t) is introduced in order to 
cover the case of a movable LOF for reasons that will be 
clear below, such as to consider the possibility to follow 
the DS movement.) 

The difference AI{t) between the intensities 2?+ and 
I?_ of the two output ports of the beam splitter, with 

P±(r,t) = i=[ii,out(r,t)±aL(r-rL(i))], (39) 



iV£;H,out(i)- 



(40) 



where the projection of the output signal field out(r, t) 
on the LOF has been introduced as the field -Eh, out (0 



(41) 



^H,out(i) -^H,out(t)+i|i,outW, 

iH.out(t) = ^ / dV al{r - TL (i))ii,out(r, i), (42) 
v N J 



with the LOF intensity denoted by 



(43) 



By using ([T]), the spectrum of the difference intensity 
fiuctuations can be written as 



Vito)^ / dTe-''"^(^H,out(i + r),£;H,out(i) 



= l+y dr e-'""^:i;H,out(t + r),SH,out(i) 

— CX3 

= l + 5out(w), (44) 

which defines the squeezing spectrum 5'out(w) of the field 
exiting the cavity. When Aout(r, t) is in a coherent state, 
V{uj) = 1 and S'out(w) = 0, which defines the standard 
quantum limit. Light is said to be squeezed at a fre- 
quency uJc when 5out('^c) < 0, and the case of complete 
noise reduction, or perfect squeezing, at ojc is signaled by 
S'out(^^c) — —1 as in this case V{lOc) = 0. 

Now the correlations of the output fields can be written 
in terms of those of the intracavity fields by using the 
input-output formalism [37| 

: il,outW,^i.out(t') :) = 27i (: A\{t),AAt') ■) (45) 



= 27i {At{t),Ai{t')) 



and then 



Sout{to) = 2j, / dr e-'^^ {6£a{t + T),d£a{t)), (46) 



where SSait) = Suit) - {£ii{t)) with 
£Hit) = AHit) + A^it), 
-4H(t)-^ / '^''^ "L(r-rL it))Ai{r,t). 



(47a) 
(47b) 
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This expression can be written in a more compact form 
as 

oo 

^out (^) = ^ / dre-'"^(5£H {t + t) 5£h (t)), (48a) 



SEnit) = (ttL {r + p{t)) I ai {r,t)), 



(48b) 



where p (t) — ri (t) — tl (t), and the change of variables 
r — > r — ri (t) has been introduced. (Remind that ri (t) 
describes the position of the dissipative structure, that 
changes because of the diffusion introduced by quantum 
fluctuations, Eqs. ((2T|) .) 

The output squeezing spectrum (|48p can now be cal- 
culated in terms of the modal correlation spectrum ([37]): 
The intracavity field fluctuations ai (r, t) can be written 
in terms of the expansion (j20p remind that Goldstone 
modes are excluded-, so that the output squeezing spec- 
trum (l48ll takes the form 



5out(^) = ^ E ("l|v,> (aL|v,) S.^iiu), (49) 
where 



("lIv,) = / d^r a£(r + p(t))i;,(r + ri (t)) + 

dV aL(r + p(t))w+(r-f ri(i)), (50) 



and N and Sij{uj) are given by Eqs. (|43|) and ([37|l . re- 
spectively. 

The projections of the eigenvectors onto the LOF 
are given by (ctLlvi), where the scalar product p3)) is 
used, and the modal correlation terms S'ij(w) correspond 
to Eq.§7^. 

Up to this point we have treated with a complete detec- 
tion of the beam, that is, we have considered an arbitrar- 
ily large detector covering the whole transverse profile of 
the outgoing field. We can wonder now on the effect of 
the detector size when it is finite, which corresponds to 
a more realistic description. 

Thus we consider now a movable detector with finite 
transverse size S, which allows to sweep the transverse 
profile of the outgoing field and study the spatial distri- 
bution of squeezing. Mathematically, the use of a finite 
size detector corresponds to limit the spatial integrations 
in Eq. to a domain S around the " detector position" 
ro, so replacing (aL|v.j), Eq. dSO]), by 



dV al {r + p{t))v, {r + n (t)) + 



{ro,S} 



dV aL(r + p(t))i;+(r + ri(t)), (51) 



where {ro,S} represents the spatial region occupied by 
the detector, and N, Eq. (gSl), by 



{ro,S} 



dV|aL(r)r 



(52) 



{ro,S} 



Finally we can compute the squeezing spectrum mea- 
sured when the finite detector is placed at Tq through 

S'out(w;ro) = — 



{ro,S} 



E ("L|vi){ro,S} ("L|vi){ro,S} Sij{uj). (53) 



As can be noted in Eq. ((53|) , the obtained level of squeez- 
ing depends both on the area and position of the detector. 



IV. INTENSITY FLUCTUATIONS SPECTRUM 

Now we apply the technique used for calculating the 
squeezing spectrum to the derivation of spectrum of the 
intensity fluctuations. For the sake of simplicity here we 
ignore the diffusive movement of the DS, i.e., we assume 
that the detector can follow such motion. 

The intensity fluctuations spectrum can be directly qb- 
served with a single photodetector and is given by |44| 



d^r d^r'x 
(5/out(r,t)5/out(r',0)), (54) 



Vi (uj) ^ dte 



where 



and the term proportional to the intensity of the outgoing 
field is 



<5/out (r , i) = /out (r , t) - ( /out (r , ) , (55) 



/out (r, t) = Al^^ (r, t) ^out (r, t) 



(56) 



By making use of ([T]) and takin g ac count of the input- 
output relation (|45|) one obtains [21'| 

Vi{uj) = Ssn[l + Si{Lo)], (57) 

where the term corresponding to the shot noise reads 



Si 



SN 



27i / dV (/i(r,t) 



and 



{ro.S} 



(58) 



(Pr dVx 
«i(r,t)(5/i(r',0) :\, (59) 
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where 



«i(r,i) =/i(r,t)- ^/i(r,i) 
/i(r,i) =il(r,<)ii(r,t). 



(60a) 
(60b) 



As stated, our interest is centered on the quantum fluc- 
tuations properties of the DS given by When such 
form of the field is considered, the expression of the in- 
tensity fluctuations spectrum, leading to first order in 
fluctuations, reads 



+ T)S£i(t)), (61) 



where 



6£i{t) = {Ai (r,t) |ai (r,i)). 



iV = / d 



\Mr)\ 



Ai (r)^(^i(r),^t(r))^ 



(62a) 
(62b) 
(62c) 



Notice that Ai (r) is the vector which corresponds to the 
classical DS, given by Eq.®. 

As can be noted, the intensity fluctuations spectrum, 
Eq. (|61|) . has the same expression as the squeezing spec- 
trum, E q.diSl) . but with the classical DS solution acting 
as LOF [2l|. So, as it occurred previously with the out- 
put squeezing spectrum, the intensity fluctuations spec- 
trum can be written in terms of the modal correlation 
spectrum ([57)) 



^^(^) = %E(^il^0(Ai|v,)5,,(c.), (63) 



where 



oo 

(Ai|v,) = J d\ [AUr)v^{r)+Ai{r)v+{r)] . (64) 



Finally, and analogously to what we did with the 
squeezing spectrum, when a finite detector of transverse 
size S positioned at Tq is considered the intensity fluc- 
tuations are given by Eq. ([63|) after replacing (Aijvi) by 
(Ai |vi)|^^ and N by ^{ro.E} so that the intensity fluc- 
tuations spectrum measured with a finite detector of size 
S placed at rg reads 

Si{uj;ro) = -J^x 



N. 



{ro,S} 



E 

1,3 



where (Ai|vi)^^^^ ^-^ and A^{ro,E} are given, respectively, 
by Eqs. ([?T|) and ([5^ when aL is replaced by Ai. 



A GENERAL RESULT ON THE SQUEEZING 
OF DISSIPATIVE STRUCTURES 



Let us assume that we can set p = in Eq. ([50)) . This 
means that the LOF can be shifted in such a way that 
it exactly follows the diffusive movement of the dissipa- 
tive structure whose squeezing is being to be measured. 
Further, let us choose the LOF — W2a;, Eq. ([TO)) , i.e., 
a LOF whose shape is Q!l — iGx- (We notice that the 
result that follows is valid for any ckl that corresponds 
to a linear combination of 'W2x and W2y.) 

By doing this it turns out that 5£yi {t) = c^x (t), see 
Eqs. (dn]) and (I48bp . Standard techniques f28| applied to 
Eq. (ED) for i^2x 



C2x 



2a; I 



(66) 



where ^2x (^) — {^2x\^) is the noise source and C2x (t) — 
(w2a;|ai), allow the computation of the stochastic corre- 
lation {SEn [t + r) SS-n (t)), that turns out to be 

{6£n{t + T)6£n{t)) =-^Nne-^^^^^l (67) 

Then, by using Eq. ([15)) we get 



(68) 



S'out i^) 



(^/27i) 



which is the main result in [25*1 . Of course the same result 
is obtained by using Eqs. ([36]), ([37]) and ([38)) . 

It is to be remarked that Eq. ([66)) is analogous to that 
derived in [i^ for the stationary phase of the hexago- 
nal mode appearing in the Kerr cavity model, and that 
it was later interpreted in as the hexagonal pattern 
transverse linear momentum. Then we can say that the 
above result means that the transverse linear momentum 
of any stationary DS appearing in the large pump detun- 
ing limit of the DOPO model is perfectly squeezed in the 
linearized theory. This is a reasonable result as it is im- 
mediately related to the fact that the transverse position 
of the DS is completely undetermined as it diffuses with 
time. 

Eq. ([55)) implies that 5 (w = 0) = — 1 what means that 
within the validity domain of the linear theory we are us- 
ing, any stationary dissipative structure sustained by the 
DOPO in the large pump detuning limit displays perfect 
squeezing at cj = when probed with the appropriate 
LOF. As Eq. ([S5)) is independent of the kind of dissi- 
pative structure and of the system parameters, it is to 
be remarked that the result is universal and independent 
of the existence of bifurcations. Let us emphasize that 
the appropriate LOF (q;l = iGx{y)) is, in principle, easily 
implementable as it corresponds to the 7r/2 phase-shifted 
gradient of the corresponding DS envelope which can be 
easily synthesized by, e.g., Fourier filtering. 

In [25| we discussed up to what extent the assumption 
p = is reasonable: It is, indeed, as the diffusion of the 
dissipative structures is very slow because of the large 
number of photons they carry, which acts as an inertial 
mass (25i] . 
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We find it important to make here a general comment 
on the mathematical technique we have presented in the 
previous sections. We must note that the linearized ap- 
proach we have presented is valid, in principle, when all 
eigenvalues are negative: In this case all fluctuations are 
damped and it is reasonable to assume that they will re- 
main small enough. Remarkably, in our case there exists 
always two null eigenvalues, which are associated with 
the Goldstone modes (see Property 1 in Subsection II 
C). Nevertheless these null eigenvalues do not make the 
linear approach invalid as the undamped fluctuations do 
not concern any particular field mode but the position 
of the dissipative structure, which is decoupled from the 
rest of fluctuations and undergoes a continuous diffusion 
(as it occurs, e.g., with the phase difference in Nu- 
merical research carried out in vectorial Kerr cavities (23| 
reinforce this confidence. In resume: In spite of having a 
null eigenvalue, we can be confident that the linearized 
description of quantum fluctuations will be reasonably 
accurate, and that a nonlinear treatment [2^ will not 
lead to dramatically different results. 



VI. SQUEEZING PROPERTIES OF THE DOPO 
BRIGHT CAVITY SOLITON 

In this second part of the article, we study in detail the 
squeezing properties of the DOPO bright cavity soliton in 
the large pump detuning limit. First, in Subsection VI. A 
we review the main properties of this solution and then, 
in Subsection VLB, we apply the theory developed in the 
first part to it. For the sake of simplicity we shall con- 
sider in this second part only the one-dimensional case, 
that is, we assume that the DOPO works embedded in a 
waveguide that avoids diffraction in the y transverse di- 
mension whilst in the x transverse dimension the system 
aperture is arbitrarily large. Moreover, we shall assume 
from now on that the LOF used in the homodyne detec- 
tion scheme can be moved, in the transverse dimension, 
in such a way that its movement exactly matches the 
diffusive motion of the bright cavity soliton, and thus we 
shall not take into account the diffusive motion of the DS 
(see for a quantitative discussion in which we show 
that this is a reasonable assumption). 



A. The DOPO bright cavity soliton 

As stated in Subsection II. A, in the limit of large pump 
detuning we are considering here, the classical description 
of the DOPO is the PDNLSE, (Eq. ^ in Appendix B). 
It is well known that the one-dimensional PDNLSE sup- 
ports two different types of localized structures (cavity 
solitons in our context): Dark cavity solitons (tanh-type 
localized structures) in the self-defocusing case a = —1, 
and bright cavity solitons (sech-type localized structures) 
in the self-focusing case cr = -fl [sl, SO, |4l|,|4i|. We shaU 
treat in the following the bright cavity soliton and thus 



we take a = +1. 

Before going on we find it convenient to review the 
main solutions of the PDNLSE. This equation has only 
two free parameters (the parametric pump parameter fj, 
and the cavity detuning Ai) as the rest of parameters 
(7]^, h, and k^) can be easily removed by normalizing 
the time and space coordinates as well as the field am- 
plitude (see Eq. (gS]) in Appendix B). When Ai < 0, 
the trivial state, Ai{x) = 0, undergoes a supercritical 
bifurcation towards a patterned state (a roll, or stripe, 
pattern) at /.t = 1. Contrarily, when Ai > the bifur- 
cation affecting the triv ial state is subcritical and occurs 
at /i = /.tQ = yT+Aj. For u, > /Iq, and positive Ai, 
dynamic patterns are found [32| . 

The bright cavity soliton has the explicit expression 

^i3s(a;) = V2/3 e*'^sech(/3a:), (69) 

with 



= Ai ± ^p? - 1, (70a) 
cos(20)=/x-\ (70b) 

This solution exists for Ai > and 1 < ^ < /Xg p3l |. 
and is stable in a wide domain of parameters, although 
for large enough Ai it becomes unstable, through a Hopf 
bifurcation, thus appearing temporal oscillations (self- 
pulsing solitons) [33|. Then the bright cavity soliton can 
undergo three different bifurcations: (i) a tangent bifur- 
cation at /.t = 1 (the BS does not exist for < 1), (ii) a 
bifurcation [i — [Iq (the BS does not exist for [i > /Iq), 
and (iii) a Hopf bifurcation for large enough Ai that 
transforms the stationary BS into a self-pulsing BS. In 
Fig. 1 we represent in the (Ai,//) plane the domains 
of existence of the different solutions we have just com- 
mented. 

It is convenient to briefly comment here on the similar- 
ities between the bright cavity soliton and the temporal 
soliton of the nonlinear Schrodinger equation (NLSE). 
The NLSE is a conservative equation that can be written 
without free parameters. Then, there is a family of sech- 
type solitons, that coexist, and for which certain quan- 
tities (such as the energy) are conserved. This is very 
different from the PDNLSE which is a nonconservative 
equation governed by two parameters. For given [i and 
Ai, the sech-type solution is unique and can be stable or 
unstable. In fact the bright cavity soliton is not a true 
soliton and its name can be misleading when compar- 
isons between these similar but actually quite different 
equations (NLSE and PDNLSE) are made. In particu- 
lar, these differences manifest in the spectra of H and C) . 

We have shown that in order to calculate the squeezing 
spectrum, given by Eq. ([49]) . one needs to evaluate the 
spectra of both £ and £t for ^^(r) = ^i,Bs(a;), 
Eq. ([55)1 . The spectrum of L in this case has been exten- 
sively studied [42l |. It consists of a continuous spectrum. 
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with eigenvalues of the form 

\s{k) = -l + .s^[^i^^^, (71a) 
Afc = Ai+sfc2, s = ±l, /je[0,oo) (71b) 

and of a discrete spectrum with eigenvalues The 
spectra of C and need to be computed numerically, 
as analytical expressions can not be derived in general, 
which we have done by adapting the Fourier method de- 
scribed in Nevertheless, some eigenvectors of the 
discrete spectrum can be computed analytically for some 
parameter sets. In particular, we have derived the four 
eigenvectors corresponding to the bifurcation occurring 
at ^ = 1. These eigenvectors are explicitly given in Ap- 
pendix B. 

A crucial point is whether the set of eigenvectors of C 
and form a (biorthonormal) basis or not. We com- 
mented above that for the similar problem of the nonlin- 
ear Schrodinger equation, describing temporal fiber soli- 
tons, the set of eigenvectors of C and do not form a 
base [ll[. As in our case analytical expressions are not 
available, our strategy has consisted in discretizing the 
transverse spatial coordinate and numerically diagonalize 
C and . Our numerical results show that the number 
of independent eigenvectors coincides with the dimension 
of the matrices, what constitutes a numerical proof that 
for this particular transverse pattern, the eigenvectors 
of C and form a (biorthonormal) basis, whereas this 
does not happen in the conservative case. We see that 
the NLSE and the PDNLSE are very different problems 
indeed. 



B. Squeezing properties of the unidimensional 
bright cavity soUton 

Now we are in conditions for studying the squeezing 
properties of the DOPO bright cavity soliton. Of course 
the amount of squeezing to be obtained will depend on 
the local oscillator field (LOF) one uses in the homodync 
measurements and, in general, on the parameter values 
(pump and detuning). We shall proceed in our study as 
follows: First, in Subsection VI.B.l, we shall particular- 
ize the general result obtained in [25] and derived above 
in Section III for the bright cavity soliton. Then, in Sub- 
section VLB. 2, we shall consider the squeezing at the bi- 
furcation points (which have been discussed in Subection 
VI. A). Finally, in Subsection VLB. 3 we shall consider the 
case of a plane-wave LOF. In this last case we discuss 
first how the amount of squeezing changes with the pa- 
rameters (for arbitrarily large detectors), and then how 
the level of squeezing depends on the size of the photode- 
tectors when these have a finite transverse dimension. 



1. Squeezing of the cavity soliton linear momentum 

In Section VI we showed that there is a special mode 
that is perfectly squeezed in the linear approximation, 
which is equivalent to saying that the transverse linear 
momentum of the DS is perfectly squeezed. We saw that 
the fluctuations of this mode are detected by choosing 
the LOF a.i^[x) — mv2x{x), that is, ai,{v) = idxAi^B,s [x), 
see Eq. (ITOl) . In Fig. 2 we represent the amplitude of this 
LOF with a dashed line together with the soliton ampli- 
tude, given by Eq. ([M)) . (As the phase factor exp (i0) 
is space independent, Eq. (|69p . it has been discarded in 
making this plot.) 

The similarity between W2x {x) and the Gauss-Hermite 
function 

Gill (x) =ie'^a;e- 5 (72) 

is immediate (in Fig. 2 the Gauss-Hermite function is 
also plotted), what suggests the use of this function, 
which is relatively easy to generate, as a LOF. In order 
to calculate the level of squeezing that would be detected 
by using the Gauss-Hermite LOF, Eq. ([7^ . we take ad- 
vantage of the general theory presented in the first part 
of this article: We expand the LOF on the basis of £t as 

= ^ (vj|q:l) Wi, (73) 

i 

excluding the Goldstone modes, and the squeezing spec- 
trum is given by 

Sontiuj) = ^ E (^»I"L> (v,|aL) Sy-(c.), (74) 

where the modal correlation matrix, Eq. (|37p , is evaluated 
from the computed spectra of £ and . The accuracy of 
the numerical method was checked by computing 5'out(i^) 
when OLi^(x) ~ W2(x), yielding an error less than 10"-'^^. 
The influence of the Gauss-Hermite LOF width and po- 
sition was already presented in [25| . and we reproduce 
it in Fig. 3 for the sake of completitude. Notice that 
the level of squeezing is quite large even when the LOF 
position or width are not well matched to those of mode 

W2(x). 



2. Squeezing at the bifurcation points 

At the bifurcations points there is, at least, one null 
eigenvalue, apart from that corresponding to the Gold- 
stone mode. This implies the existence of a mode, dif- 
ferent to W2(a;), whose eigenvalue reaches its minimum 
value. This mode is expected to be perfectly squeezed in 
the linear approximation. The squeezing of these modes 
is the equivalent to the usual squeezing at the bifurca- 
tion points that has been repeatedly studied in a large 
number of nonlinear cavities 13, m| . 
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We have seen that the bright cavity soUton can un- 
dergo three different bifurcations at /i = 1, /i = /ip, and 
/i — /ijjBj see Fig. 1. Let us consider first the bifur- 
cation occurring at = 1 (the analytic expression of 
the mode that has A = — 2 is given in the Appendix 
B, mode W3(x)). Then, by taking ol\,(x) ~ W3(a;), one 
obtains the squeezing spectra that we have represented 
in Fig. 4 for Ai — 1.2 and three values of /i. No- 
tice that 5out {cl! — 0) ~ —1 ior fj, ~ 1 (full line); for 
/i = 1.0001 (dashed line) the maximum squeezing does 
not reach —1. Furthermore, as fi departs from unity the 
maximum squeezing does not occur at w = 0, but at a 
slightly different frequency. For fj, = 1.01, the maximum 
level of squeezing is close to —0.75. That is, the squeezing 
degrades quickly as the system departs from this bifur- 
cation. A similar behavior is exhibited at the bifurcation 
occurring at /i = /ip = \/l + and we shall not enter 
into quantitative details. 

Let us finally consider the Hopf bifurcation. For 
/i — /ijjB there are two eigenmodes, let us denote them 
as WHB+ and whb-, for which the eigenvalues read 
X± = —2 zbiijjHB (whb is the oscillation frequency at the 
Hopf bifurcation), i.e., there are two most damped eigen- 
modes which are expected to exhibit maximum squeez- 
ing. But in the Hopf bifurcation there cannot be a LOF 
that matches these squeezed modes: Remind that the 
LOF vector was defined as q:l = (Q^L,Qf£)"^, and at 
the Hopf bifurcation the eigenvectors with Re A = —2, 

WHB± = {wHB±,w^B±) ' do not verify Whb± = "^HBi- 
This can be appreciated in Fig. 5 where we have rep- 
resented the imaginary parts of w^^^ and w^q_^ (also 
the real parts, not shown, are quite different). This fact 
makes it impossible to find a LOF such that (q^lIvrb) — 
1. Then perfect noise reduction is never achieved at the 
Hopf bifurcation. 

The amount of squeezing attainable at the Hopf bi- 
furcation is represented in Fig. 6 where the squeezing 
spectrum has been represented for two choices of the 
LOF: ai^{x) = WHB+(a:), in dashed line, and ai,{x) = 
WHB+(a^) +whb-(2;), in full line. Interestingly the max- 
imum squeezing is larger for the latter choice. Notice 
also that, in contrast to the other bifurcations and as 
it is well known, at the Hopf bifurcation the maximum 
level of squeezing is reached at a frequency different from 
u! = 0, which corresponds to the self-pulsing frequency of 
the bright cavity soliton. 

3. Squeezing with a plane-wave local oscillator 

We consider now the particular case of the squeezing 
properties of the bright cavity soliton when measured 
in a homodyning experiment, probed with a plane- wave 
LOF, which is of obvious relevant interest as this is the 
simplest LOF. 

At a first stage we deal, as in the preceding subsec- 
tions, with a complete beam detection, that is, the de- 
tector is assumed to completely cover the transverse ex- 



tension of the outgoing quantum cavity soliton. Fig. 7 
shows the squeezing spectra obtained for five different 
values of the pump fi and Ai = 1.2. The maximum of 
squeezing is achieved at frequency lo — irrespective of 
the pump value (which is obviously limited to 1 < /i < 
= ^1 + A^, see Fig. 1). A high degree of squeez- 
ing, sustained for a large range of parameters setting, is 
obtained. This is more clearly seen in the inset of Fig. 
7, where S{uj — 0) is depicted as a function of pump for 
Ai = 1.2. 

We have also calculated the influence of the detun- 
ing: In Fig. 8 we plot the same as in Fig. 7 but for 
fixed fi = 1.2 and as a function of detuning Ai. Again 
large squeezing levels are obtained in all the domain of 
existence of the cavity soliton (for detunings larger than 
those represented in Fig. 8, the cavity soliton becomes 
Hopf unstable as already commented). 

We focus now on the homodyning detection when us- 
ing finite size detectors. When a detector, with size S 
and positioned at x — xq is used, the squeezing spec- 
trum detected is given by Eq. (|53p . particularized to the 
BS solution. In the particular case of a plane-wave LOF, 
the calculations are considerably simplified, as is con- 
stant. We note, as well, that the phase of the plane-wave 
LOF is a free parameter, which is allowed to vary in order 
to obtain the maximum level of squeezing. 

In Fig. 9, we represent the spatial distribution of 
squeezing, for fi — Ai = 1.2 and ai lo — 0, when the 
finite size detector is displaced across the transverse di- 
mension. We consider three different values of the de- 
tector size E = Ax/ (3 normalized to the BS width f3, 
see Eq. (|70ap . as indicated in the figure. In the three 
cases the phase of the plane-wave LOF (not shown) has 
been chosen in order to obtain the maximum squeezing at 
each position of the detector in the transverse plane. A 
clear conclusion can be extracted from these plots: The 
smaller is the detector size, the smaller is the squeezing 
level one attains. This is to be expected as the smaller 
the photodetector is, the larger influence of high spa- 
tial frequency modes, i.e., with the smaller photodetector 
the influence of all sort of nonsqueezed modes is larger 
than with the larger photodctectors which filter out high 
frequency modes. One more conclusion one can extract 
from the figure is that the squeezing level is larger, for all 
detector sizes, at the center of the soliton {x = 0), i.e., 
the soliton is more squeezed than the vacuum (for large 
X the soliton amplitude tends to zero and the squeezing 
is due to the squeezed vaccum). 

Finally we focus on the influence of the detector size 
on the degree of squeezing reached when the detector is 
placed at the exact center of the soliton. Fig. 10 shows 
the maximum squeezing reached in these conditions at 
frequency uj — when a detector of normalized size E 
is considered. Several features can be appreciated in the 
figure. As a general trend, the squeezing level degrades 
as the detector size is reduced, i.e., the maximum level of 
squeezing is found for a very large detector. But, inter- 
estingly, for 1.5 < E < 3, the squeezing decreases with 
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the increase of the detector size. Obviously, for detec- 
tors whose size falls in this region what is happening is 
that the fluctuations being detected come from both the 
BS and the trivial solution (below threshold emission) 
existing far from the detector. In fact we see that the 
phase of the local oscillator that optimizes the squeez- 
ing level (represented with a dashed line, left vertical 
axis), which is almost insensitive to the detector size for 
S < 1, changes rapidly in this detector width region. 
Then, for large detectors, S > 10, the monotonic increase 
of squeezing with the detector width is recovered and the 
optimum phase for the LOF is again almost insensitive 
to the detector size. 



VII. CONCLUSIONS 

In this article we have developed a general theory for 
the analysis of linearized quantum fluctuations of optical 
dissipative structures generated in wide aperture nonlin- 
ear optical cavities. Although we have done this explic- 
itly for the special case of the degenerate optical para- 
metric oscillator in the large pump detuning limit, our 
method can be easily generalized to any other nonlinear 
cavity. The method consists, in short, in expanding the 
fluctuations in the biorthonormal base that the eigenvec- 
tors of the linear deterministic operator of the linearized 
Langevin equations of the system. This technique al- 
lows, in particular, the identification of a special mode 
which is perfectly squeezed (in the linear approximation). 
The perfect squeezing occurs irrespective of the nonlin- 
ear cavity parameter values, and the special mode can 
be identified with the transverse linear momentum of the 
dissipative structure that is being emitted. It must be 
emphasized that the existence of this squeezed mode is a 
genuine transverse effect, i.e., it is associated to the sym- 
metry breaking introduced by the existence of dissipative 
structures. 

Then we have applied this theory to the study of 
squeezing of a particular dissipative structure, namely, 
the bright cavity soliton. In particular we have analyzed 
the squeezing occurring at the different bifurcations that 
this dissipative structure can undergo at the same time 
that the appropriate LOF for detecting it in each situa- 
tion. Then we have studied the squeezing detected when 
a plane wave LOF is used and analyzed its dependence 
on the parameter values. Finally we have considered also 
finite size detectors. We have shown that for large detec- 
tors the squeezing level is large almost independently of 
the system parameters. For finite size detectors, we have 
analyzed the spatial distribution of squeezing. 
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VIII. APPENDIX A 

In this Appendix the DOPO model used along the pa- 
per is derived. The Hamiltonian describing the DOPO 
in the interaction picture is given by [2^ 



(75) 



where 



H,,,, ^nJ2^^fd\ Ai (A„ ~ llV^) i„, (76) 



governs the free evolution of the intracavity fields in the 
paraxial approximation. 



i^int = I / d^^ 



io(il) -4 [A, 
describes the nonlinear interaction, and 



H, 



h hrri 



(77) 



(78) 



accounts for the coherent driving. In the above expres- 
sions In = c/ ^'Zujn'^n the diffraction length for the 
field An, Aq = (wo - 2cjs) and Ai = (loi - w^) /-f-^ 
are the (adimensional) pump and signal detuning param- 
eters, respectively, V = d"^ jdx^ -t- jdy^ is the trans- 
verse Laplacian operator, and g is the (real) nonlinear 
coupling coefficient, given by 



5 



2(27rn)3 V £oi. ' 



(79) 



with ~xP'^ the relevant nonlinear susceptibility matrix el- 
ement, n the common value of the refractive index of the 
crystal at pump and signal wavelengths (a type I DOPO 
is considered), and the thickness of the crystal along 
the resonator axis. 

From the above Hamiltonian one obtains the master 
equation governing the evolution of the density matrix p 
of the intracavity modes, 



d 



1 



dt^ ih 
where the Liouvillian term 



H,p 



Ap, 



(80) 



Ap= J2 [ 'i''^ [2i„(r,t)piUr,t) 

^ — n 1 



n=0.1 



pAi{r,t)An{r,t) - An{r,t)Ai{r,t)p\ , (81) 

models the coupling between the system and the external 
reservoir through the output mirror. 

Passing to the the generalized P representation one 
can transform the master equation (j80p into an equiva- 
lent Fokker-Planck equation for a quasiprobability den- 
sity (denoted by P), following standard methods (see, 
e.g. [201), the result being: 
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dAt 



(-l.LlAt + ^-At' £l) + 

-71^1^1 - gAtAo) + 

-1^LIA+ - gAiA+) + 
92 



_9_ 

d 
dAt 
g ( 9^ 



-4n 



^+ ] F(A). (82) 



2 VSy^P*^ dAX^- 
In the above expression A = [Aa, A^ , A\, , and 



2v72 



(83) 



In its turn, a Fokker-Planck equation, here Eq. ([82ll . 
can be transformed into an equivalent classical-looking 
set of stochastic differential equations (so called Langevin 
equations) via Ito rules |28|] . In our case they read 





= IqLqAo 






(84a) 


^A^ 


= 70^0-4^ 




U f* 
r '--in J 


(84b) 




= 7iii^i 




+ V5A?7(r,t), 


(84c) 


^A^ 

at' 




+ gAiA^ 


+ y.g<77+(r,t). 


(84d) 



with 77 (r, t) and ?y^(r, t) independent, real white Gaussian 
noises of zero average and correlations given by 

(77+(r',t'),77(r,i)) =0, (85a) 

(ry+(r,t),77+(r',i')) = (5(r - r')(5(i - t'), (85b) 
{rj{r,t),7jir',t'))=S{r-r')6{t-t'). (85c) 

Note that, if Af is interpreted as A* and the noise 
terms are ignored (classical limit), Eqs. coincide 
with the classical equations for a planar DOPO [29j . 

Equations ([51)1 are already set in a convenient way to 
apply the method we develop in this paper. However we 
shall use a simpler form of them obtained in the limit of 
large pump detuning, defined by |Ao| >> 
which allows the adiabatic elimination of the pump field 
faEm as 



Ao = Ao 



ad 



7o^o 
7o^o 



A+ - A+ - ^ (f* ^A+^ 



(86a) 
(86b) 



We note that these expressions correct some typos ap- 
pearing in the corresponding expressions in ^25]. The 
remaining Langevin equations for the signal field are ob- 
tained by substitution of the limit solution (|86p into 



By taking, without loss of generality, the external pump 
field amplitude as a purely imaginary quantity whose 
imaginary part has the same sign as the pump detun- 
ing Aq, i.e. = icr \£in\ with cr — sign Aq, they become 
Eqs. ([2]) in Sec. Ill Al which is the model we shall con- 
sider. 



IX. APPENDIX B 

In the classical limit {A^ — > Al, and 77,77^ 0), Eqs. 
([2|) reduce to 

^ A (r, t) = 7i (ii^i +fiAl+i^ \Ai I' a) , (87) 

Li = —(1 -t- iAi) + ilfV^, which is a version of the so 
called parametrically driven, nonlinear Schrodinger equa- 
tion (PDNLSE), which is a universal model for pattern 
formation in parametrically driven systems (see, e.g., (35j 
for a list of systems described by the PDNLSE) . A con- 
venient normalization of that equation is obtained by in- 
troducing the following dimensionless quantities: time 
T = jit, spatial coordinates {X,Y) = (x/li,y/li), and 
field i/j ~ Ai/k, so that Eq. ([57|) becomes 



drip = ^l4*-i^+iAl)^P+i {d\ + 4) ijj+ia ^. 



For the bright cavity soliton, Eq. ((69|) . operators C 
and take the form 



where 



^-^ Al CI y ^ - \Ai Ci )^ 



2(3 



(89a) 



£1 = -(l + iAi) + iV^ + q — ) sech^(/3a;), (90a) 



A - M + i ( ^) e^"^ sech^ [px) 



and 



= Al ± vV^, 
cos(20) — fi~^. 



(90b) 



(91a) 
(91b) 



Then, for = 1 one can calculate the discrete eigen- 
values analytically. The result is that the discrete eigen- 
vectors of C and with null eigenvalue are 



Vl 



Wi 



-5r 



-s 



e 



{x + iT) 

[x ~ iT) e-"l' J ' 

[p^ +i{xT~ 1)] e*"^ 



W4 = P^^^S 



(92a) 
(92b) 
(92c) 
(92d) 
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and that the discrete eigenvectors of £ and C"^ with eigen- In the above expressions Vi is the Goldstone mode, Eq. 

value A = —2 are ([7|), and we have introduced the quantities 

^2 = -^P' S^_^^^_^^^^^.,, j, (93a) 

-1/2 c-r/" ro.M . fP 



W2^~iP 'ST\^_^_i^j, (93b) 5 = Y'^secM/Jx), T = /3tanh(/3a;). (94) 

V3 = ^^3S ) , (93c) 

wc, = -B^^S ( y\ ^ * ^^'^ ^ '^'^} ^ f93d) Let us finally notice that for /i = 1, it turns out that 

V [/3 /32 = Ai andcos(2<^) = l. 



[1] P. D. Dmmmond and Z. Ficek (eds.), Quantum Squeezing 

(Springer, 2004). 
[2] R. Loudon and P. L. Knight, J. Mod. Opt. 34, 709 (1987). 
[3] P. Meystre and D. F. Walls (eds.), Nonclassical Effects 

in Quantum Optics (American Institute of Physics, New 

York, 1991). 

[4] N. Korolkova, G. Leuchs, R. Loudon, T. C. Ralph, and 
Ch. Silberhorn, Phys. Rev. A 65, 052306 (2002). 

[5] V. V. Dodonov and V. I. Man'ko (eds.). Theory of Non- 
classical States of Light, (Taylor and Francis, London, 
2003) 

[6] D. F. Walls and G. J. Milburn, Quantum Optics, 

(Springer Verlag, Berlin, 1994). 
[7] S. Reynaud, C. Fabre, and E. Giacobino, J. Opt. Soc. 

Am. B 4, 1520 (1987). 
[8] M. D. Reid and P. D. Drummond, Phys. Rev. Lett. 60, 

2371 (1988). 

[9] P. D. Drummond and S. J. Carter, J. Opt. Soc. Am. 4, 
1565 (1987). 

[10] H. A. Haus and Y. Lai, J. Opt. Soc. Am. B 7, 386 (1990). 
[11] V. V. Kozlov, IEEE J. Sel. Topics Quant. Electron. 9, 
1468 (2003). 

[12] N. Treps and C. Fabre, Phys. Rev. A 62, 033816 (2000). 
[13] E. Lantz, T. Sylvestre, H. Maillotte, N. Treps, and C. 

Fabre , J. Opt. B: Quantum Semiclass. Opt. 6, S295 

(2004). 

[14] E. M. Nagasako, R. W. Boyd, and G. S. Agarwal, Opt. 

Express 3, 171 (1998). 
[15] A. Mecozzi and P. Kumar, Quantum Semiclass. Opt. 10, 

21 (1998). 

[16] L. A. Lugiato, M. Brambilla, and A. Gatti, Optical Pat- 
tern Formation, Adv. Atom. Mol. Opt. Phys. 40, 229 
(1999). 

[17] L. A. Lugiato, A. Gatti, and E. Brambilla, J. Opt. B: 
Quantum Semiclass. Opt. 4, S176 (2002). 

[18] M. I. Kolobov, Rev. Mod. Phys. 71, 1539 (1999). 

[19] L. A. Lugiato and G. Grynberg, Europhys. Lett. 29, 675 
(1995). 

[20] L. A. Lugiato and A. Gatti, Phys. Rev. Lett. 70, 3868 
(1993). 

[21] A. Gatti and L. A. Lugiato, Phys. Rev. A 52, 1675 
(1995). 

[22] P. D. Drummond and K. Dechoum, Phys. Rev. Lett. 95, 
083601 (2005). 

[23] R. Zambrini, M. Hoyuelos, A. Gatti, P. Colet, L. A. 
Lugiato, and M. San Miguel, Phys. Rev. A 62, 063801 



(2002). 

[24] D. Gomila and P. Colet, Phys. Rev. E 66, 046223 (2002). 

[25] I. Perez- Arjona, E. Roldan, and G. J. de Valcarcel, Eu- 
rophys. Lett. 74, 247 (2006) 

[26] A. Gatti, H. Wiedemann, L. A. Lugiato, I. Marzoli, G. L. 
Oppo, and S. M. Barnett, Phys. Rev. A 56, 877 (1997). 

[27] P. D. Drummond and C. W. Gardiner, J. Phys. A 13, 
2353 (1980). 

[28] C. W. Gardiner and P. ZoUer, Quantum Noise (Springer, 
2000). 

[29] G. L. Oppo, M. Brambilla, and L. A. Lugiato, Phys. Rev. 

A 49, 2028 (1994). 
[30] S. Longhi, Phys. Scr. 56, 611 (1997). 
[31] S. Trillo, M. Haelterman, and C. Sheppard, Opt. Lett. 

22, 970 (1997). 

[32] S. Longhi and A. Geraci, Appl. Phys. Lett. 67, 3060 
(1995). 

[33] M. Bondila, I. V. Barashenkov, and M. M. Bogdan, Phys- 

ica D 87, 314 (1995). 
[34] I. V. Barashenkov, N. V. Alexeeva, and E. V. 

Zemlyanaya, Phys. Rev. Lett. 89, 104101 (2002). 
[35] G. J. de Valcarcel, I. Perez- Arjona, and E. Roldan, Phys. 

Rev. Lett. 89, 164101 (2002). 
[36] M. J. Collet and D. F Walls, Phys. Rev. A 32, 2887 

(1985). 

[37] M. J. Collet and C. W. Gardiner, Phys. Rev. A 30, 1386 
(1984). 

[38] N. Alexeeva, I. V. Barashenkov, and G. P. Tsironis, Phys. 

Rev. Lett. 84, 3053 (2000). 
[39] S. Fauve and O. Thual, Phys. Rev. Lett. 64, 282 (1990). 
[40] J. W. Miles, J. Fluid Mech. 148, 451 (1984). 
[41] C. Elphick and E. Meron, Phys. Rev. A 40, 3226 (1989). 
[42] I. V. Barashenkov, M. M. Bogdan, and V. I. Korobov, 

Europhys. Lett. 15, 113 (1991). 
[43] N. V. Alexeeva, I. V. Barashenkov and D. E. Pelinovsky, 

Nonlinearity 12, 103 (1999). 
[44] R. J. Glauber, Phys. Rev. 130, 2529 (1963); ibtd. 131, 

2766 (1963). 

[45] A. Gatti and S. Mancini, Phys. Rev. A 65, 013816 (2001). 



15 



FIGURE CAPTIONS 

Fig.l. Schematic representation of the different pat- 
terns one finds in the PDNLSE, Eq. (gHl). The hne fi^B 
has been ploted after the numerical results of [33l] . 

Fig. 2. Amplitude of the bright cavity soliton, Gauss- 
Hermite mode (GHi), and the mode W2 (in dashed line). 

Fig. 3. Squeezing level (at the labeled frequencies) 
displayed by the bright cavity soliton when a Gauss- 
Hermite mode (GHi) or mode W2 (see text) is used as a 
LOP. In (a) a GHi of width ^ is used. In (b) the LOP is 
displaced x from the CS center. Parameters are Ai = 1, 
and fi — 1.2. 

Fig. 4. Squeezing spectra close to the tangent bifur- 
cation at /i = 1.0. The detuning parameter is Ai = 1.2 
and the pump values are indicated in the figure. 

Fig. 5. Imaginary parts of w^^^ (full line) and w^^^ 
(dashed line). See text. 

Fig. 6. Squeezing spectrum at Hopf bifurcation. The 
dashed line corresponds to using one of the two most 
damped modes (with Re A = 0) at this bifurcation as 
LOP. The full line corresponds to using as LOP the sum 



of these two modes. 

Fig. 7. Squeezing spectra obtained when using a plane- 
wave LOP. The detuning parameter is Ai = 1.2 and the 
pump value is indicated in the figure. In the inset the 
maximum squeezing at w = is represented as a function 
of pump for the same detuning value. 

Fig. 8. Squeezing spectra obtained when using a plane- 
wave LOP. The pump parameter is /i = 1.2 and the de- 
tuning value is indicated in the figure. In the inset the 
maximum squeezing at w = is represented as a function 
of detuning for the same pump parameter value. 

Fig. 9. Spatial distribution of the squeezing when 
probed with a plane-wave LOP and measured with finite 
size detectors of normalized size S = Ax/ (3 (the values 
are marked in the figure) . The phase of the LOP has been 
chosen in order to obtain the maximum noise redution at 
each detector position. Parameters are /i = Ai = 1.2. 

Fig. 10. Maximum squeezing Sout a.t cu — 0, obtained 
with a finite size detector centered at the BS center, as a 
function of the detector size S = Ax/p. Parameters are 

= Ai = 1.2. 
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